NOVEMBER  1976. 


SUFFIELD  MEMORANDUM  29/76 


NT  1 5 REPRODUCTION 

,oto>m5Sj»  ii  J**1ass",h> 


PROJECT  NO.  14-60-06 


DfPBNCK  RISC  ARCH  ISTABUSHMINT  SUNIILD  J RALSTON  S ALMRTA 


COMPUTER  SIMULATION  OF  ANOMALIES  CREATED  BY  AN  ARRAY  OF 
FERROUS  OBJECTS  IN  THE  EARTH'S  MAGNETIC  FIELD  (U) 


R.U.  Hemrlch 


UNCLASSIFIED 


UNCLASSIFIED 


TABLE  OF  CONTENTS 


Table  of  Contents 


List  of  Figures  11 

Abstract  ill 


Introduction 


Theory 


Program  Operation 


A.  Inputs  and  Operation 


B.  Output 


Appendix  10 


(1) 

UNCLASSIFIED 


' — ' WIWUI  mi  WIW'i, . i jmi  jg  mi 

UNCLASSIFIED 

LIST  OF  FIGURES  Following 

Page 

1.  Typical  Anomaly  due  to  a Ferrous  Object  1 

2.  Anomaly  Components  3 

3.  Earth's  Field  Vector  4 

4.  Plan  view  of  Array  4 

5.  Anomaly  From  an  Array  of  Five  Objects  6 

6.  Anomaly  at  a Field  Inclination  of  45°  6 

7.  Typical  Anomaly  from  Two  Objects  in  South-East  Alberta  6 


(11) 

UNCLASSIFIED 


UNCLASSIFIED 


DEFENCE  RESEARCH  ESTABLISHMENT  SUFFIELD 
RALSTON  ALBERTA 


SUFFIELD  MEMORANDUM  NO.  29/76 


COMPUTER  SIMULATION  OF  ANOMALIES  CREATED  BY  AN  ARRAY  OF 
FERROUS  OBJECTS  IN  THE  EARTH'S  MAGNETIC  FIELD  (U) 


R.W.  Heinrich 


ABSTRACT 

^ This  report  describes  a FORTRAN  IV  computer  simulation  of  a three- 
dimensional  array  of  ferrous  objects  and  Its  effect  on  the  earth's 
magnetic  field.  The  program  output  Is  a plot  of  the  magnetic  anomaly 
amplitude  versus  distance  along  a path  similar  to  that  recorded  by  a 
total  field  magnetometer  passing  over  the  array.  The  array  size,  object 
mass,  direction  of  the  earth's  field  and  the  altitude  and  direction  of 
the  path  of  the  magnetometer  are  controlled  by  the  user.  One  hundred 
and  one  computer  words  of  storage  are  required  for  each  pass. 

This  study  was  done  as  part  of  a Range  Clearance  Technology  task 
leading  to  the  design  of  equipment  for  the  detection  of  unexploded  ordnance 
found  on  military  ranges. 
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FERROUS  OBJECTS  IN  THE  EARTH'S  MAGNETIC  FIELD  (U) 

by 

R.W.  Henri ch 

INTRODUCTION 

Previous  work  Involving  magnetic  anomalies  due  to  a dipole  has  been 
restricted  to  the  signal  obtained  from  a single  object  (1).  The  purpose 
of  this  study  was  to  derive  a general  expression  for  the  effect  of  a 
ferrous  object  on  the  earth's  magnetic  field  and  adapt  It  for  computer 
use  to  allow  the  calculation  of  the  anomaly  due  to  a three-dimensional 
array  of  these  objects.  This  work  would  not  only  allow  the  simulation 
of  the  unexploded  ordnance  scattered  about  a drop-zone,  but  would  permit 
the  modeling  of  shrapnel  lying  near  the  surface,  or  the  modeling  of  an 
area  of  Iron-rich  soil  of  high  susceptibility. 

Figure  1 shows  a plot  of  anomaly  magnitude  due  to  a single  ferrous 
object  as  seen  by  a total  field  magnetometer  scanning  over  the  object. 
Where  an  array  of  these  objects  occurs.  It  Is  difficult  to  visualize  the 
total  effect  that  they  will  have  on  the  earth's  field. 
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In  order  to  create  a smooth  plot  of  the  total  effect  due  to  an 
array  of  objects  It  Is  necessary  to  determine  the  effect  at  a large  number 
of  points  on  the  path.  In  this  simulation,  101  points  were  used  on  a 10" 
10"  plot.  Where  the  array  contains  100  objects  the  result  Is  10,100 
calculations  per  pass.  For  this  reason  the  problem  lends  Itself  well  to 
a computer  solution. 

THEORY 

An  object  In  the  earth's  magnetic  field  may  cause  a distortion  of 
the  field  depending  on  the  permanent  and  Induced  magnetization  of  the 
object.  It  Is  this  distortion,  or  anomaly,  which  indicates  the  presence 
of  an  object. 

Permanent  magnetization  is  determined  by  the  history  as  well  as 
composition  of  the  object  and  can  add  to  or  cancel  the  effect  of  the 
Induced  component.  For  this  study  only  induced  magnetization  is  con- 
sidered. 

Induced  magnetization  is  a function  of  the  Intensity  of  the  ambient 
field  and  the  magnetic  susceptibility  of  the  material  in  the  object: 

I = KH 

Where  I * Induced  magnetization  per  unit  volume 
K * magnetic  susceptibility 
H * earth's  magnetic  field 
The  magnetic  moment  of  the  object  Is 
M * IV 

Where  V * volume 
Therefore  M « KHV 

The  magnetic  moment  for  typical,  spherical,  man-made  iron  objects  is 

5 6 

10  to  10  cgs  units  per  metric  ton  (1). 
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The  anomaly  created  by  an 
the  distance  to  the  object  and 
T = CM 


object  at  a point  Is  a function  of 
Is  of  the  form 


Where  C * constant 

n * constant 

r * distance 

Objects  which  have  high  length  to  diameter  ratios  tend  to  behave  as 

monopoles,  especially  If  magnetized  along  their  length.  That  Is,  their 

field  lines  radiate  straight  out  from  their  ends  and  the  resulting  anomaly 

2 

Is  Inversely  proportional  to  r , or  n ■ 2.  In  reality,  they  are  dipoles 
which  have  their  poles  far  apart  (1). 

In  this  study,  objects  of  Interest  have  low  1/d  rat'  and  behave 
as  dipoles.  In  this  case  n Is  equal  to  3,  and  the  anomaly  Is  Inversely 
proportional  to  r3. 

A magnetic  anomaly  due  to  a dipole  Is  sometimes  expressed  as  two 
components.  One  component  Is  due  to  the  field  off  the  end  of  a dipole.  In 
this  case  constant  C Is  equal  to  2.  The  other  component  Is  from  the  side  of  the 
dipole  where  C Is  1,  as  shown  In  reference  2.  The  first  component  is 
projected  onto  the  line  joining  the  object  of  detection  and  the  point  of 
Interest  and  becomes  the  radial  component  Tr.  The  other  Is  projected  onto  the 
perpendicular  to  It  and  becomes  the  tangential  component  T*  as  shown  In 
Figure  2. 
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Direction  of  the 
earth's  field. 


Figure  2 Anomaly  Components 
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The  two  components  are  (1) 
Tr  = 2M  cos  <* 
r3 

T«  = -M  sin  « 


A three-dimensional  co-ordinate  system  is  used  to  determine  the  angle 
which  is  the  angle  between  the  earth's  field  vector  and  the  line 
between  object  and  magnetometer.  The  vector  describing  this  line  is 
projected  onto  a unit  vector  representing  the  earth's  field.  The  angle  <* 
can  then  be  determined  from  this  projection.  The  earth's  field  vector  VI 
is  defined  here  (Figure  3). 

A1  = cos  A cos  B 

B1  = sin  A cos  B 

Cl  = sin  B 

where  VI  = 1 

The  vector  V2  describing  the  direction  and  distance  from  the 
magnetometer  to  any  object  in  the  array  is  defined  below  (Figure  4). 

A2  = 1X2  - y sin  0 + x cos  0 
B2  = IY2  + y cos  0 + x sin  0 
C2  = H 

where  H is  the  height  of  the  magnetometer  above  the  top  2-dimensional 
array  of  objects. 

and  V2  * V A22  + B22  + C22 
Vectors  VI  and  V2  can  be  described  as 
VI  » Ali  + Blj  + Cl k 
V2  * A21  + B2j  + C2k 
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The  dot  product  Is  used  to  project  one  vector  on  another: 
7T«  V2  * V1V2  cos  « 
or 

cos  « * V1*V2 
VI V2 

but 

V1*V2  ■ A1A2  + B1B2  ♦ C1C2 
and 

V1V2  » Y A22  + B22  ♦ C22 
Therefore 

cos  « » A1A2  ♦ B1B2  + C1C2 
VA22  ♦ B22  ♦ C22 

and 

sin  « -Vi  - cos2  « 

We  can  now  find  Tr  and  T«  since 

Tr  - 2Hcos  « 
r3 

and 


T«  ■ -Ms In  « 
r3 

These  components  are  then  projected  Into  the  earth's  field  vector.  This 
Is  done  because  a total  field  magnetometer  measures  the  component  of  an 
anomaly  very  nearly  In  the  direction  of  the  earth's  field  (1).  The 
direction  of  the  earth's  field  Is  easily  determined  and  remains  quite 
constant  for  a given  area. 
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Tz  ■ Tr  cos  « + T«  sin  « 

Tz  is  the  anomaly  detected  by  the  magnetometer  for  one  object. 

The  computer  program  adds  the  anomalies  for  each  object  In  the 
array  and  produces  plots  like  those  shown  in  figures  5-7.  Figure  5 
Is  a west- to- east  scan  over  an  array  of  objects  5x5x1.  Each  object 
appears  In  a scan  this  close  to  the  array  whereas  a scan  such  as  that 
shown  In  Figure  6 Is  far  enough  above  the  array  of  2 x 2 x 3 to  hide  the 
detail  of  each  object.  Figure  6 Is  also  an  example  of  a south-to-north 
scan  In  a location  where  the  Inclination  of  the  earth's  field  Is  45°.  All 
Positive  Inclinations  are  below  the  horizontal  In  this  program.  Figure  7 
shows  a scan  where  the  Inclination  Is  73°,  approximately  the  angle  In 
south-eastern  Alberta.  The  symmetry  which  occurs  In  a vertical  field 
is  noticeably  absent  at  this  angle.  A magnetic  moment  of  5 x 10s  cgs 
units/metric  ton  has  been  used  to  produce  these  curves. 

PROGRAM  OPERATION 

A.  Inputs  and  Operation 

The  program  Inputs  and  Its  basic  operation  are  explained  here. 

Refer  to  the  main  program  listing  In  the  Appendix  for  Input  formats. 

PROBS  number  of  Independent  plots  done  in  a run. 

THETA  angle  of  rotation  of  the  array  In  degrees  about 

a vertical  line  through  the  reference  point  (Figure  4). 

A declination  of  the  earth's  field  In  degrees  when  the 

magnetometer  heading  Is  north  (Figure  3). 

B Inclination  or  dip  of  the  earth's  field  In  degrees. 
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Figure  5 Anomaly  From  An  Array  of  Five  Objects 


Figure  6 Anomaly  at  a Field  Inclination  of  45 
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XP  horizontal  distance  In  centimeters  from  the  center 

of  the  array  to  the  nearest  point  of  approach  on 
the  magnetometer  path. 

S horizontal  distance  In  centimeters  from  the  center 

of  the  array  to  the  start  of  the  magnetometer  path. 

F horizontal  distance  In  centimeters  from  the  center 

of  the  array  to  the  end  of  the  magnetometer  path. 

N01  number  of  objects  In  the  array  In  the  X-dl recti on 

(Figure  4). 

N02  number  of  objects  In  the  array  In  the  Y-dlrectlon 

(Figure  4). 

H vertical  distance  In  centimeters  from  the  top  two- 

dimensional  array  to  the  nearest  point  on  the 
magnetometer  path. 

M distance  In  centimeters  between  objects  In  the  array 

in  the  X-dl recti on. 

N distance  In  centimeters  between  objects  In  the  array 

In  the  Y-dlrectlon. 

WT  weight  of  the  objects  In  the  array  In  kilograms. 


NUM  maximum  number  of  objects  In  the  array  In  the  vertical 

direction. 

P vertical  distance  In  centimeters  between  each  two- 

dimensional  horizontal  array. 

The  main  program  reads  In  the  data  cards  In  the  order  given  above. 
Then  the  reference  points  (XXR,  YYR)  are  defined  relative  to  the  center 
of  the  array. 

Subroutine  SCAN  Is  called  to  do  the  anomaly  calculations.  These 
calculations  are  all  made  relative  to  the  reference  point.  The  distance 
1X2  from  the  reference  point  to  the  closest  point  on  the  magnetometer 
, path  Is  printed  out  prior  to  SCAN.  The  subroutine  sets  the  location 
of  the  simulated  magnetometer  according  to  the  user  Instructions.  Then 
at  101  equally  spaced  points  along  this  path  It  determines  the  anomaly. 
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The  anomaly  at  a point  Is  determined  by  finding  the  angle  « 
and  the  distance  to  each  object  and  using  the  expression  derived  earlier 
to  obtain  the  total  anomaly  from  the  array.  This  data  Is  stored  In  computer 
core  In  a two-dimensional  array  called  AA  which  Is  dimensioned  to  AA  (101,4). 
It  has  room  for  four  magnetometer  passes  on  one  graph.  Three  are  shown  In 
Figure  6.  The  first  line  plotted  Is  for  the  top  two-dimensional  array.  A 
new  line  Is  then  plotted  as  an  additional  two-dimensional  array  Is  added 
below  the  last  and  the  effect  of  all  present  arrays  Is  calculated.  In 

Figure  6 three  lines  are  shown  Indicating  a three-dimensional  array,  three 
layers  vertically,  the  last  line  being  calculated  for  the  entire  array. 

Array  size  Is  set  by  the  user.  As  the  data  Is  being  stored,  the  maximum 
anomaly  among  the  303  points  Is  determined  for  use  later.  When  storage 
Is  complete,  control  Is  returned  to  the  main  program. 

Subroutine  SKALE  Is  called  and  uses  the  maximum  anomaly  point  found 
In  SCAN  to  set  an  appropriate  scale  for  the  y-axls  of  the  plotter  output. 

Again  control  Is  returned  to  the  main  program.  The  line  printer 
will  output  the  maximum  anomaly  point  T and  the  maximum  scale  height 
XMAX  which  were  found  In  SCAN  and  SKALE. 

The  grid  Is  set  with  four  tick  marks  per  inch  for  ten  Inches  on 
both  the  X and  Y axis  of  the  plotter. 

Subroutine  ROUND  Is  called  to  round  off  the  X axis  scale  If  necessary. 
The  main  program  then  Instructs  the  plotter  to  print  the  X axis  numbers, 
one  per  Inch  as  required  by  the  scale,  and  the  X axis  title  Is  printed. 

The  same  procedure  occurs  for  the  Y axis. 
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The  rest  of  the  main  program  plots  the  data  stored  In  array  AA. 

When  the  plot  is  finished  the  plotter  pen  moves  to  a new  location 
ready  to  begin  another  graph  If  Instruction  PROB  requires  It. 

B.  Output 

The  main  program  output  Is  a graph  of  the  magnetic  anomaly  due  to  an 
array  of  ferrous  objects.  At  the  side  of  the  graph  Is  a listing  of  the 
Input  parameters  to  the  program.  The  anomaly  scale  In  gammas  Is  printed 
along  the  vertical  axis.  The  distances  In  centimeters  printed  on  the 
horizontal  axis  give  the  horizontal  component  of  the  distance  from  the 
geometric  center  of  the  object  array  to  the  simulated  location  of  the 
magnetometer.  The  program  Input  lists  the  desired  start  and  finish 
distance  from  the  magnetometer  to  the  array  center.  These  values  are 
transformed  Into  the  horizontal  scale  of  the  graph. 

As  mentioned  in  SECTION  A,  the  magnetometer  heading  is  always 
assumed  to  be  north.  However,  by  adding  to  the  declination  angle,  the 
effective  heading  can  be  changed.  Also,  the  entire  array  may  be  rotated 
about  Its  reference  point  (SECTION  A)  to  compensate  for  this  new  heading 
(Figure  4). 

The  distance  XP  Indicates  the  closest  distance  that  the  magnetometer 
comes  to  the  array  center.  When  XP  Is  zero,  the  detector  Is  passing 
directly  over  the  geometric  center  of  the  array. 
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RFAL  M,N 

DIFFUSION  AA (101*4) 

RFAO  I 7 » 7) PROPS 

7 FORMAT (13) 
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RFAO  ( 7 • 4 ) TmFT A 

4 FORMAT ( F9.1 ) 

THFTA  • THFTA  * 2.#  3.14/960. 

RFADI2*  6 1 A*B 
6 F0RMATIFS.1.FS.1I 
A ■ A ♦ 90. 

A • A • 2. #3. 14/960. 

P • " • 2. *9. 14/960. 

RFA0(2*  8IXP 

8 FORMAT U 10 1 
READ(2.  10 )S*F 

10  FORmAT(|10*I10) 

RFAO  <2*12 )N01 *N02 
12  FORMAT  119*19) 

RFAO  (2.14JH 
14  FORMAT  ( F10.2 I 
2 m M 

RFAO  (2*16)M,N 
16  *ORMAT  |F10.2*F10,2) 

RFAO  (2*18 )WT 
IP  FORMAT  (P10.2) 

C CONVERT  to  MFT9IC  TON 5 

WT  ■ WT/100C .0 
®FAO  (2*20) MUM 
20  FORMAT  It 31 
RFAO  I2.22)P 
77  FORMAT  (F10. 21 

XXR  • m# ( NOl-1 ) /?•  ♦ M/2, 

YVR  • N*(N02-l)/2.  ♦ M/2. 

TL  • SORT  (XXP**2  ♦ YYR**2) 

AN  • ATAN  ( YVR /XXR ) 

ANTm  ■ AN  ♦ THFTA 
XU*  COS  (ANTH)  * -L 
CALL  ROWNO  (XU*ll 


PAGF  P 


PAGF  9 

PO  • PO  ♦ ,1  * xmax 

OK.-XMAX+FlOAT ( J J ) / 10.0FXMAX 
11  CONTI NUf 

YX  ■ -FLOAT ( S*F)/10.-FLOAT(L> 

XXX. ,89  »XMAX-XMAX 

CALL  FCMAR  ( YX.  XXX ,C. 1 *0. 1 • 1 .9* I 
WRITF  (7,34) 

34  FORMAT  ( • ANOMALY  ( GAMMA I * I 

M • 2 

WRITF  19*  31 
3 FORMAT  COM 
00  13  X*1,NUM 

LL-1 

MM.-FLOAT(L) 

29  A1*0. 

00  19  Il*l*X 

19  A1*A2  +AA ( LL  *1 1 ) 

IF(Al-XMAX)  17*19,21 
17  IFCAlTXMAX)  29*19,19 
23  A1*-XMAX 
GO  TO  19 
21  Al*  KMAX 
GO  TO  19 

19  CALL  FPL0Tl-2.HM.Al) 

LL*LL*1 

MM.HM4-FL0AT<  ( 5*F)/100) 

I F ( MM— ( I V— 2*L ) I 29*29*27 
27  CALL  FPLOT  (♦l.MH.Al) 

13  CONTI NUF 

A « A/2. /3. 14*360, 

N • B/2./3.14*36C. 

CALL  FCMAR  (<-L  ♦ ll./X)*  0.94  • XMAX * 0.1*0. 1*0.0) 
WRITF ( 7*  36)M 

96  FORMAT  ( 'HEIGHT *M'11X'. SF9. 2. 1X*CM. • ) 

CALL  FCMAR  ( ( — L ♦ ll./X).  0.48  • XMAX*  0.1*0. 1*0.0) 
A * A - 90. 

WRITF  ( 7*38 ) A 

38  FORMAT  C OFCL I NAT  I ON.A • 6X • * ' *F9. 1 *2X • DEG. M 

CALL  FCMAR  ( I — L ♦ ll./X).  0.4 2 * XMAX,  0.1*0. 1*0.0) 
WRITF (7.  40 ) P 

40  FORMAT  C JNCIINAT!0N.P'6X'»'.F8.2*2X'DFG,M 

CALL  FCMAR  ( ( -L  ♦ ll./X),  0.36  • XMAX*  0.1*0,l*0.0) 
WRITF  (7.42IXP 

42  FORMAT  ('DISTANCE, XP*PX»** ,I6*4X'CM.M 

CALL  FCMAR  ((-L  ♦ ll./X),  0.30  * XMAX,  0.1*0. 1.0,0 
THFTA  • THFTA  / 2./  3.14*360. 


1X2  • XP  - I F 1 X ( XU ) 

WRITF  13.1) 

1 FORMAT  Cl'  I 

WRITF(3*  24)1X2 
24  FORMAT ( ' ','1X2  *'*I10) 

RAY  ■ S!N(ANTH)*Tl 
CALL  ROWNO  (RAY * 1 ) 

L ■ S - I F I X ( RAY ) 

IV  • S ♦ F ♦ L 

CALL  SCAN  (M,N*M.N01.N02.WT.AA.THFTA.A*B.IX2*L.NUM.Z.P.T.S*F) 
CALL  $XAWE  (T.XMAX) 

WRITF  (3.26)T,XMAX 

26  FORMAT ( • O'  * 'T  * • .F19.6 *3X , • XMAX  *'.F13.6) 

Y*1 , /XMAX* 10, /2, 

X*10./FL0AT(S*F) 

PF*-FLOAT(L) 

PU*— XMAX 

CALL  SCALF(X.Y,PF.PU) 

XXY..029*FLOAT(S*F| 

CALL  FGRI0(C*P?.PU*XXv,40) 

CALL  FGRID  (0,PF.r.0.XXY,<*0) 

VY  • XMAX/40. 

CALL  FGRID(1.PF.PU.YY.80) 

PO*-.39*FLOAT(  5*p ) /10.— F_OAT(L ) 

XO  • -S 

00  R If  • 1*11 

VYY«— ,04*XMAX— XMAX 

CALL  FCMAR  (PO*  YYY  *.1*.1*0.) 

WRITF  (7,281X0 
28  FORMAT (19  I 

PO*PO*FLOAT(  S*F)/1C. 

XO*XO*(  S*F)/10 
9 CONTINUE 

XY..44*FL0AT(  S*F)-FLOAT(Ll 
XX*— .OP* XMAX— XMAX 

CALL  FCMAR  I XY,  XX  *0, 1 *0* l *0.0 ) 

WRITF  (7*30) 

30  FORMAT ( 'POSITION  (CM.)') 

PO*-.OOS*XMAX-XMAX 

OX  ■ -XMAX 

00  11  JJ  • 1*21 

XVX.-FLOAT(  S*F)*. 07-FLOAT (L) 

CALL  FCMAR  ( XYX *P0* *1 *• 1 *0*0 ) 

CALL  ROWNO  (OX .2) 

WRITE  (7,92>OK 
32  FORMAT ( Fp. 1 ) 
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WR I TF ( 7*  44IM 

44  FORMAT  CX-OTR.  5PACIN6*M'9X'*»  *F7.0*3X'CM.' I 

CALL  FCMAR  ( ( -L  ♦ ll./X)*  0.24  * XMAX*  0.1*0. 1*0.0) 
WRI TF ( 7,46 ) N 

46  FORMAT  ( 'V— DIR.  SPACING.N* 3X ••• *F7,C*3X'CM. • ) 

CALL  FCMAR  ((-L  ♦ ll./X),  0.18  * XMAX*  O.l.C. 1*0.0) 
WRITE  (7,5P)  P 

98  FORMAT  ('VERTICAL  SPACI \G,P • *1X • • • ,FQ.2 * IX . • CM. • ) 
CALL  FCMAR  ((-L  ♦ ll./X),  0.12  * XMAX*  0.1,0. 1*0.0 
WRITF  (7.48IN01 

48  FORMAT  ( * ARRAY  W I DTH *NOl • 4X • ■ ' * I 6 ) 

CALL  FCMAR  ((-L  * ll./X)*  0.C6  • XMAX*  0.1*0.1*0.0) 
WRITE  (7,90)  N02 

90  FORMAT  ('ARRAY  LENGTH *N02 • 3X '• ' *16) 

CALL  FCMAR  M-L  * ll./X),  0.00  * XMAX*  0.1*0.i*0.9) 
WRITE  (7,36)NUM 

96  FORMAT  ('ARRAY  DEPTH *NUM» ,4X '■ ' * 16 ) 

CALL  FCHAR  ( C -L  ♦ ll./X), -0.06  • XMAX,  0.1*0,l*0.0) 
WRITE  (7*92)  THFTA 

92  FORMAT  ( »R0TATI0N*THFTA'5X'*' *F7.e*3X'0EG. ' ) 

CALL  FCMAR  ((-L  ♦ ll./X), -0.12  * XMAX,  0.1*0. 1*0.0) 
C CONVERT  TO  XG. 

WT  • WT  • 1000.0 
CALL  ROWNO  ( WT*3 ) 

WRITF  (7, 34 )WT 

94  FORMAT  ( •WFIGMT.WT'lOX1**  *F9.2*lX'XG.' ) 

APR  ■ 1.7*1  FLOAT ( S*F ) )— FLOAT (L ) 

CALL  FPLOT  ( *3 • APP •— XMAX ) 

7 CONTINUE 
CALL  FXIT 
END 
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